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Vortex structure in density wave scenario of pseudogap 

Maciej M. Maska and Marcin MierzejewskQ 
Department of Theoretical Physics, Institute of Physics, University of Silesia, 40-007 Katowice, Poland 

We investigate the vortex structure assuming the ci-density wave scenario of the pseudogap. We 
discuss the profiles of the order parameters in the vicinity of the vortex, effective vortex charge and 
the local density of states. We find a pronounced modification of these quantities when compared 
to a purely superconducting case. Results have been obtained for a clean system as well as in 
the presence of a nonmagnetic impurity. We show that the competition between superconductivity 
and the density wave may explain some experimental data recently obtained for high-temperature 
superconductors. In particular, we show that the rf-density wave scenario explains the asymmetry 
of the gap observed in the vicinity of the vortex core. 
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I. INTRODUCTION 

Recent scanning-tunneling-microscopy (STM) exper- 
iments give a new insight into the electronic structure 
of the vortex cores in the high-temperature supercon- 
ductors (HTSC), and, more generally, into the problem 
of the interplay between superconductivity and a mag- 
netic field. In type-II superconductors, the flux pene- 
trating the vortex suppresses the superconducting order 
parameter locally in the core. It was shown in 1964 
by Caroli at al.^ that there should be bound states 
around the vortex core in an isotropic s-wave supercon- 
ductor. Hess was the first who has experimentally con- 
firmed the existence of core states or core excitations in 
NbSe2 superconductor 0. These states manifest them- 
selves in the local density of states (LDOS) as a large 
peak near the Fermi surface. Similar experiments car- 
ried out for HTSC have shown a rich structure of the 
vortex states. Namely, Maggio-Aprile at al. have found 
a splitting of the central peak in YBa2Cu3 07_x|3]. The 
resulting peaks are separated by about 11 meV. On the 
other hand the numerical solution of the Bogoliubov-de 
Gcnnes (BdG) equations for c?-wave superconductivity 
(DSC) indicates the absence of such a splitting. 0| It was 
later explained by Franz and Tesanovic that there are no 
truly bound states in d-wave superconductors and all the 
states are extended with continuous energy spectrum(5i]. 
They tried to explain this discrepancy by introducing an 
additional dxy component of the superconducting order 
parameter. The resulting dx-2-y2 -\- idxy gap is nodeless 
and leads to an exponential decay of the core states. Gen- 
erally, the presence of the bound states is expected for 
any nodeless gap. Although, dx^-y^ is commonly ac- 
cepted as a dominant symmetry of the order parameter 
an additional component, that closes the nodes in the 
X — ±y directions, can be induced by the external maj 
netic field y. However, results presented in Ref. 
suggest that qualitative agreement with the STM exper- 
iments can be obtained also for a pure d-wave super- 
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conductivity, provided that the magnetic field is strong 
enough. Moreover, the presence of the additional compo- 
nent may not lead to an essential modification of LDOS. 
The presence of the core states in BSCCO, unlike YBCO, 
is not well established. Small peaks at ±7 meV have been 
reported in Ref. S|. On the other hand, the core states 
have not been observed in other STM experiments |^. 
Instead, a gaplike structure has been found in the center 
of the vortex core. It has been identified as a pseudogap 
that evolves smoothly into the superconducting gap away 
from the vortex core. As the origin of the pseudogap is 
still under debate, investigation of the vortex structure 
can shed a new light on this problem. 

The pseudogap has been investigated with the help of 
various experimental techniques like: angle-resolved pho- 
toemission [T l[ [l ^ . intrinsic tunneling spectroscopy 
111 0: NMR [lllil, infrared ^3 and transport 
[l8j measurements. Aside from other approaches a d- 
density-wave (DDW) state has recently been proposed 
to explain the pseudogap HI IH [p. The DDW sce- 
nario has been investigated to verify whether it actually 
applies to the pseudogap regime. For the details we re- 
fer to considerations concerning the transport properties 
13 , nonmagnetic _2^ and magnetic J23 impu- 
rities as well as the phonon self-energy [2^ |29| | in the 
DDW phase. According to this hypothesis the pseudo- 
gap opens due to condensation of electron-hole pairs with 
a nonzero angular momentum (^ — 2). As a result there 
are staggered fluxes originating from the orbital currents, 
which alter from one plaquette to the neighboring one. 
The DDW order breaks the time reversal, rotational and 
translational invariance, preserving combination of arbi- 
trary two of them. There exists a similar approach that, 
however, does not break the translational symmetry [22j |. 

As the DDW and DSC orders compete (1^ |3^ , one 
may expect an enhancement (or appearance) of the DDW 
gap near the vortex core, where the superconductivity is 
suppressed. It is possible due to the insensitivity of the 
DDW order to the magnetic fleld jsj. In the present 
paper we investigate the vortex structure assuming the 
DDW scenario of the pseudogap. Our approach explains 
the tunneling spectra obtained in the vicinity of the vor- 
tex core in BSCCO0. We analyze the vortex structure in 
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a wide range of doping level and find crucial differences 
between under- and overdoped regimes. It was shown 
in Ref. that depending on the carrier concentration 
one may expect pure DDW, mixed DDW+DSC, and pure 
DSC ordering. Our results suggest the possibility of an 
additional phase, where the DDW order occurs only in- 
side the vortex cores. In the framework of the SU(2) 
slave-boson theory this possibility has been investigated 
by Kishine et al. |33] . As the vortices are often pinned at 
impurities , we investigate the electronic structure also 
in such a case. Finally, we discuss a possible coupling be- 
tween the staggered currents and the antiferromagnetic 
order, that has recently been observed in the vicinity of 
vortices. Isa 



II. THE MODEL 



where 



W,, = (A) 4c,. - r;,, (A) c],c..) (7) 

is the DDW amplitude. Note, that in contradistinction 
to the superconducting order parameter magnetic field 
exphcitly enters Wij j31l] . 

The mean-field Hamiltonian can be diagonalized with 
the help of the Bogoliubov-de Gennes (BdG) equations. 
Namely, we introduce a set of new fermionic operators 
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Our starting point is an effective Hamiltonian that de- 
scribes a system with coexisting DDW and DSC orders 

M 

H = iJkin + ^^DSC + ^^DDW- (1) 

The kinetic part of the Hamiltonian is given by 

where t is the nearest-neighbor hopping integral and /i is 
the chemical potential, cj^ {cia) creates (annihilates) an 
electron with spin cr at site i. According to the Peierls 
substitution [35| the magnetic field enters the kinetic part 
of the Hamiltonian through a phase factor r\ij (A) : in the 
presence of magnetic field the hopping between sites i and 
j is accompanied by acquiring of an additional phase, 
given by 



where 



Vjl 



Uil 
Vll 



Here, the single particle Hamiltonian is given by 



(8) 



~tS,+s.j + i-iyW,,] 7y„- (A) - My , (9) 



and both the order parameters are determined self- 
consistently by: 



(-1) 



i^DDW 



^Im [uiiu*ir],j (A) 



V 2 



+ vuv*,rj*^ (A)] tanh (^^^ , (10) 
Vbsc J2 (^Uiiv*i + v*iUji) tanh (^^) -(H) 



,, (A)=exp {jj^^^-d^ 



(3) 



Husc is the nearest-neighbor pairing responsible for the 
d-wave superconductivity 



where 



Ay = ^(caCjT 



(4) 



(5) 



is the superconducting order parameter. As we do not 
specify the mechanism responsible for pairing, Vdsc is 
assumed to be field-independent . |3^ The DDW state oc- 
curs due to 

HuDW = J2 i-'^yW^i (A) clc,., (6) 



In the presence of magnetic field the phase factor 77 ex- 
plicitly appears in the expression for the DDW order pa- 
rameter [Ea. HlU|) ]. It differs the BdG equations used in 
this paper from those obtained in Ref. |25l |. 

In order to compare the numerical results with STM 
data we have calculated the local density of states 
(LDOS) 



{e) = -J2 IV' [E^ + \vu\'f {Ei + e)] 



(12) 

where /'(e) is the derivative of the Fermi distribution 
function /(e) = [exp(e/fcr) -|- 1] . Pi{e) is proportional 
to the local differential tunneling conductance that could 
be measured in STM experiments. We have also calcu- 
lated the local electron concentration 

^2j2Wa\^fiEi) + \va\^fi~Ei), (13) 



that determines, e.g., the effective charge of the vortex. 
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III. RESULTS 

We have carried out calculations for 35 x 35 square 
lattice with one superconducting flux quantum piercing 
this area. We have taken the nearest-neighbor hopping 
integral as the energy unit and assumed Vddw = 1-6, 
Vdsc = 1-4- In the absence of magnetic field the phase di- 
agram of HTSC can qualitatively be reproduced for such 
values of the interaction strengths (we refer to Ref-psf 
for the details). 

A. Structure of vortex 

We have analyzed the vortex structure for different val- 
ues of the occupation number. For small concentration 
of holes, 6 < 0.04 ((5 = 1 — n), the DDW order dominates 
and superconductivity is completely suppressed. There- 
fore, there is no vortex for such an occupation number 
and the DDW state is homogeneous. For larger doping 
the DDW and DSC orders coexist. For a doping slightly 
higher than S — 0.04 both the phases coexist in the whole 
sample, as depicted in Fig. ^ The presented magnitudes 
of the order parameter are defined by 

= i {A.u+i + A,,,_x - A, - A»,,_y) (14) 

DDW and DSC orderings compete with each other [l^ . 
Therefore, reduction of one of them enhances the other 
one. As the DDW order is hardly affected by the mag- 
netic field it becomes strongly enhanced in the vortex 
core where the superconductivity is suppressed. In or- 
der to prove that this mechanism is responsible for the 
enhancement of the DD W order in the vortex core we 
have plotted < $ >= y'($pDW^2 ^ ($psc)2 ^^^^ ^.^^ 

EJ. Note, that this quantity is almost constant. 

Farther increase of doping {S > 0.1) destroys the DDW 
order in a homogeneous svstem. j25|| However, this order- 
ing can be restored in the vortex core. Such a case is pre- 
sented in Fig. O The local coexistence of DDW and DSC 
takes place for the doping up to (5 ~ 0.16. Therefore, due 
to the suppression of superconductivity, magnetic field 
significantly enhances the doping regime where the DDW 
order occurs. For higher doping the DDW order is not 
restored in the vortex core. It originates from the fact 
that for sufficiently strong doping DDW docs not exist 
even in the absence of superconductivity (Vbsc — 0)- 

Another interesting feature, that originates from the 
presence of the DDW order is related to charge of the vor- 
tex. For small doping the vortex charge is negative (i.e., 
the electron concentration in the vortex core is higher 
than outside the vortex, see Fig. ^). As the doping 
increases the electron concentration in the vortex core 
decreases and for 6 > 0.16 the vortex core becomes pos- 
itively charged. The inversion of the vortex charge oc- 
curs when the DDW order disappears. It originates from 




FIG. 1: The vortex structure obtained for /i = —0.22 {S = 
0.075) and kT = 0.05. The upper panel shows the DSC order 
parameter, whereas the lower one corresponds to the DDW 
order. 



the changes of LDOS: the DDW gap opens in the cen- 
ter of the band and enhances the density of states below 
the Fermi level. This inversion may also be induced by 
changing of temperature, provided that the DSC transi- 
tion temperature is larger than those of the bare DDW 
state (overdoped region). Namely, at low temperature 
the DDW order is restored at the vortex core and, there- 
fore, vortex is negatively charged (see Fig. 2t>)- At suffi- 
ciently high temperature the DDW order does not occur 
and the vortex is positively charged (Fig. ^). Exper- 
imental observation of this feature would certainly sup- 
port the DDW scenario of the pseudogap. 
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FIG. 2: < $ ^($pDW-)^ _^ ($°^°)^ calculated for the 
same parameters as in Fig. 



B. Local density of states 

The first STM measurements of the vortex core in 
BSCCO did reveal neither the existence of the core states 
nor the zero bias conductance peak (ZBCP) expected for 
the pure DSC state . The presence of the core states in 
BSCCO has later been reported in Ref. Q. The discrep- 
ancy between both the experiments remains, however, 
unexplained. In our approach the DDW and DSC order 
parameters have nodes in the same parts of the Fermi 
surface and the bound states do not occur (inclusion of 
components with other symmetries is needed). However, 
we show that the DDW scenario explains the qualitative 
features of LDOS reported in Ref. 0. In Fig. |S1 we 
present the LDOS for different distances from the vortex 
center. For such occupation number both the DDW and 
DSC gaps are present also away from the vortex. When 
approaching the vortex center the DSC gap, that exists 
at the Fermi level, vanishes and the DDW gap becomes 
strongly enhanced. The DDW gap opens in the center 
of the band and, therefore, is responsible for the asym- 
metry of the LDOS inside the vortex. In particular, the 
peak at positive bias shifts outwards when approaching 
the center of vortex, whereas the peak at negative bias 
does not move. Such an asymmetry has recently been 
observed in BSCCO '^fj . For larger doping this feature is 
even more pronounced since there is no DDW gap away 
from the vortex and the Fermi level is much below the 
center of the band. 

The LDOS at the vortex center evolves smoothly with 
doping. This is depicted in Fig. |^ The lowest curve cor- 
responds to the case when DDW and DSC orders coexist 
in the bulk material. The following four curves are ob- 




FIG. 3: The same as in Fig. 1, but for fi = -0.38 {5 = 0.16). 



taincd for DDW and DSC orders coexisting in the vortex, 
whereas the two topmost curves are obtained for a purely 
DSC vortex. For small doping the DDW gap does not 
allow for the formation of the ZBCP. The DDW gap de- 
creases with increasing doping and the ZBCP gradually 
develops in the overdoped regime. 

In the overdoped regime a similar evolution of LDOS, 
connected with the vanishing of the DDW gap, can also 
be caused by the increase of temperature. Fig. \7\ shows 
such a situation for S — 0.16. At low temperature the 
DDW gap strongly suppresses the ZBCP, however a rem- 
nant of this peak is visible at the Fermi level. At tem- 
perature A:T = 0.1 the DDW gap closes and ZBCP be- 
comes strongly enhanced. Farther increase of temper- 
ature destroys superconductivity, ZBCP disappears and 
the normal-state van Hove singularity appears in the cen- 
ter of the band (the uppermost curve, kT = 0.25). 
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FIG. 4: The electron concentrations in the vicinity of vortex. 
The upper panel (a) has been obtained for {5 = 0.1 and kT = 
0.05. The middle (b) and lower (c) panels show the results 
for kT = 0.05 and kT = 0.1, respectively. In cases (b) and 
(c) S = 0.16. 




FIG. 5: LDOS calculated for 5 = 0.075 and kT = 0.05 for var- 
ious distances from the vortex center. The uppermost curve 
corresponds to the vortex center, whereas the lowest one to 
the distance of 10 lattice constants from the center. 




-1.5 



FIG. 6; Doping dependence of LDOS calculated at the 
vortex center at kT = 0.05. The curves from the 
bottom to the top correspond to the dopings: S = 
0.075,0.1,0.12,0.15,0.16,0.17,0.18. The dashed line shows 
the Fermi level, whereas the arrows indicate the center of 
band. 



Decay of the DDW order, when moving away from the 
vortex core, is doping dependent. In the underdoped case 
the DDW order parameter decreases monotonically from 
the maximum in the vortex core to the value of the bulk 
DDW+DSC state. In the overdoped case the DDW or- 
der vanishes outside the vortex, however the decay in not 
monotonic but rather oscillating. There exist lines where 
the DDW order changes sign, what corresponds to re- 
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FIG. 7; Temperature dependence of LDOS calculated at 
the vortex center for 5 — 0.16. The curves from the bot- 
tom to the top correspond to the temperatures: kT — 
0.01, 0.05, 0.10, 0.15, 0.20, 0.25. The broadening has been as- 
sumed to be temperature independent. 



verted circulation of the staggered currents (see Fig. |SJ). 
The resulting pattern posses a fourfold symmetry and has 
a checkerboard-like modulation. However, the amplitude 
decays rapidly and disappears at the distance of several 
lattice constants. The size of the checkerboard squares 
is comparable to the coherence length (a few lattice con- 
stants) and increases with decreasing doping. Similar 
pattern has recently been observed in the STM imaging 
of slightly overdoped BSCCO The staggered flux 

alone may be to weak to explain the antiferromagnetic 
order observed in the vicinity of the vortex core |33| . 
However, presence of DDW order may stabilize the an- 
tiferromagnetic ordering of spins. Here, two mechanisms 
can be taken into account: (i) a coupling between the 
staggered magnetic field and spins, {ii) changing of the 
local electron concentration towards half-filling, where 
the antiferromagnetism is stable (see Fig. This 
tempting hypothesis should, however, be verified within 
a separate study (e.g., one can extend the approach pre- 
sented in Ref. |37|). This problem is currently under 
investigation. 



C. Vortex pinned at impurity 

As a large number of vortices is pinned at impurities 
we have also analyzed the vortex structure and LDOS 
in such a case. Investigation of the LDOS in the vicinity 
of a nonmagnetic impurity has recently been proposed 
as an experimental test of the DDW scenario pa 13 ■ 
Since the DDW order is strongly enhanced in the vortex 
core one can expect a nontrivial modification of the vor- 



FIG. 8: The contour plot of the DDW order parameter in the 
vicinity of the vortex calculated for S = 0.16 and kT = 0.01. 



tex structure. In order to account for the presence of a 
nonmagnetic impurity we have extended the Hamiltonian 



[Eq. ©] by 



imp 



(16) 



where U is the strength of the impurity located at Rq. 
Due to the presence of the impurity the single particle 
Hamiltonian in the BdG equations acquires an additional 
term 



(17) 



Fig. ini shows the DDW order parameter around the 
vortex that is pinned at a weak impurity with U = 2. 
The DSC order parameter is absent in the vortex core 
and, therefore, weakly affected by the impurity. Con- 
trary to this, the DDW order is is completely suppressed 
at the impurity site, provided that U is large enough. 
Then, the DDW currents are absent on the bonds, which 
link the impurity site with its nearest neighbors. As a 
result the DDW order is suppressed also at the neighbor- 
ing sites and has a maximum on a ring surrounding the 
vortex center. The spatial structure of the pinned vortex 
is reflected in the LDOS, (see Fig. I1U|) . At the impu- 
rity site the most of the spectral weight is transfered to 
larger energies. At the nearest neighbor site the DDW 
order is still suppressed by impurity, whereas DSC gap 
is negligible as we are in the vortex core. This allows for 
the formation of the ZBCP, as can be inferred from the 
second curve in Fig. ^| With increasing distance from 
the impurity both the order parameters develop and the 
ZBCP disappears. Here, one can observe the asymmet- 
ric gap, as discussed before for a clean superconductor. 
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FIG. 9: DDW order in the vicinity of vortex that is pinned 
at impurity. We have used 5 = 0.1, kT=0.05 and (7 = 2. The 
vortex center is at the left corner of figure. 

Then, the DDW order parameter decays and the pure 
DSC order parameter achieves its bulk value (the upper- 
most curve in Fig. IIUI) . 

The above results show important qualitative differ- 
ences between the structure of vortices uninfluenced and 
influenced by impurities. Assuming the DDW origin of 
the pseudogap, the ZBCP should occur only in the latter 
case. Therefore, experimental verification of the DDW 
scenario requires the ability to independently map the 
vortex and impurity locations. Such an experimental 
technique has recently been described in Ref. 

IV. SUMMARY AND REMARKS 

In conclusion, we have investigated the the vortex 
structure in the system with coexisting DSC and DDW 
correlations. Changing the carrier concentration we have 
found four phases: (i) pure DDW phase, {ii) DDW and 
DSC orders coexisting in the whole system, (m) DSC 
order existing in the bulk and DDW order only in the 
vortex core, and (iv) pure superconducting vortex. The 
system goes from the phase (i) to (iv) when the doping 
increases. The actual positions of the phase boundaries 
are model-dependent, and may change when the elec- 
tron correlations are taken into account. In the case (ii) 
the DDW order parameter is strongly enhanced in the 
vicinity of the vortex core. It is due to the competition 
between DSC and DDW orderings. The presence of the 
DDW gap in the vortex [phases {ii) and {Hi)] does not 
allow for the formation of ZBCP. In this case, the evolu- 
tion of the DSC gap into DDW one, when approaching 



the vortex center, provides a natural explanation of the 
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FIG. 10: LDOS calculated at various distances from the vor- 
tex center. The model parameters are the same as used in Fig 
9. The lowest curve corresponds to the vortex center, whereas 
the uppermost one to the distance of 5 lattice constants from 
the center. 

gap asymmetry reported in Ref. This asymmetry 

originates from the fact that the DDW and DSC gaps 
open at different energies. 

As the DDW order suppresses the ZBCP, this peak 
can be observed in the phase {iv), where the system is 
purely superconducting. However, one can expect the 
occurrence of the ZBCP when the vortex is pinned at 
impurity. This is due to the absence of the staggered 
currents in the close neighborhood of the impurity. 

As DDW and DSC gapes have nodes in the same di- 
rections, the core bound states do not occur in our ap- 
proach. However, an additional component of the order 
parameter, for which the gap is nodeless, may lead to the 
splitting of the zero-bias peak in LDOS. 

The presented modification of the vortex structure 
originates predominantly from the competition between 
DSC and DDW orders. A direct coupling of the DDW or- 
der to the magnetic field is of minor importance. There- 
fore, one can argue that the presented results remain 
valid also for other non-superconducting order param- 
eter, that competes with DSC and is weakly affected by 
the magnetic field. 
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